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CONVERSION  OF  A  CYLINDRICALLY  CONFINED  SURFACE  EXPLOSIOls 
INTO  A  ONE-DIMENSIONAL  BLAST  WAVE 

I.  INTRODUCTION 

In  connection  with  the  Dense  fock  (Close-Spaced  Basing)  scheme,  it  has 
been  pointed  out  by  Latter  that  surface  environments  following  a  multiburst 
attack  could  be  considerably  more  difficult  to  survive  than  those  resulting 
from  single  bursts.  Detonation  simultaneously  over  all  the  shelters  of  a 
closely  spaced  array  surrounds  the  vicinity  of  an  individual  shelter  with  high 
pressures,  allowing  the  explosion  to  vent  only  in  the  upward  direction. 

Hence,  at  least  in  the  interior  of  the  array,  extremely  high  pressures  are 
maintained  for  longer  times,  with  the  result  that  total  impulses  applied  to 
ground  structures  are  nuch  greater  than  for  a  solitary  explosion. 

Latter  proposed  considering  the  following  idealized  situation:  identical 
explosions  are  initiated  simultaneously  in  an  infinite  regular  hexagonal  array 
at  (or  near)  a  uniform  level  plane  surface  (Fig.  1).  Construct  a  vertical 
plane  transverse  to  the  line  connecting  every  pair  of  neighboring  centers,  at 
the  midpoint  of  the  line.  By  symmetry,  the  six  planes  surrounding  an 
explosion  represent  perfectly  reflecting  surfaces  and  together  form  a 
hexagonal  parallelepiped  (a  cylinder  having  hexagonal  cross  section)  with  a 
vertical  axis.  This  in  turn  can  be  approximated  as  a  right  circular  cylinder 
with  the  same  axis  and  cross  sectional  area.  A  two-dimensional  r-z  code  can 
be  used  to  solve  this  idealized  problem  for  surface  and  air  bursts.  It  is 
clear,  however,  that  asymptotically  the  flow  becomes  one-dimensional,  evolving 
into  a  shock  tube  solution.  In  the  limit  where  this  asymptotic  state  has  been 
Manuacript  approved  November  9,  1982. 
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reached  and  where  many  weapon  masses  of  air  have  been  swept  over,  the  one¬ 
dimensional  form  of  the  Sedov  point  blast  solution  should  apply.  The 
pressure  on  the  ground  as  a  function  of  time  t  is  given  by 


,  W  2/3 

p  =  k  p  — -  , 

^  o ' p  At  ' 
o 


(1) 


where  pQ  is  the  air  density,  W  is  the  weapon  yield,  A  is  the  cylinder  area, 

and  k  is  the  function  of  the  adiabatic  index  y  which  is  <  1.  This  is  to  be 
compared  with  the  corresponding  result  for  the  case  of  a  spherical  free-field 
expansion,  which  is  of  the  form 


p  -  **p0  (7T3) 

0 


2/5 


(2) 


As  may  be  seen,  the  pressure  given  by  Eq.  (l)  decreases  in  time  more 
slowly  than  that  of  Eq.  (2).  Although  the  idea  behind  this  so-called  "bomb- 
in-a-can"  model  is  simple  and  transparent,  it  leaves  several  important 
questions  unanswered.  First,  how  long  (how  many  shock  reflections)  does  it 
take  to  reach  the  asymptotic  state,  i.e.  ,  when  does  Eq.  (l)  become  a  good 
approximation?  Second,  to  what  extent  is  the  prediction  (l)  dependent  on 
idealizations  of  the  model,  e.g.  ,  perfect  symmetry,  perfect  simultaneity, 
neglect  of  scouring  and  en*- raining  of  dust,  etc?  Finally,  what  are  the  long¬ 
term  tactical  consequences  in  this  model  of  a  multiburst  scenario — 
specifically,  what  can  be  said  about  the  shape  of  the  plume  produced,  how  mich 
dust  is  raised  up  and  where  does  it  go,  what  sort  of  turbulence  is 
established  and  what  are  the  consequences  for  radar  transmission,  etc? 

In  this  report  we  describe  a  calculation  which  was  carried  out  to 
investigate  the  first  question,  regarding  the  approach  to  the  one -dimensional 
solution.  We  have  used  a  version  of  the  DNA  one  kiloton  (l-kton)  standard  to 
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initialize  a  surface  burst  equivalent  to  18  scaled  megatons  (Mton)  in  a 
cylinder  of  radius  900  feet.  We  find  that  the  flow  is  close  to  one¬ 
dimensional  after  ~20  ms  (the  time  required  for  the  blast  to  reflect  from  the 
cylindrical  wall  and  return  to  the  origin,  considerably  less  than  might  have 
been  expected. 

In  the  following  section  we  describe  in  detail  the  calculations  and  the 
results  obtained,  and  conclude  with  a  brief  discussion  of  what  we  intend  to  do 
in  follow-on  work. 


2.  TWO  DIMENSIONAL  CALCULATION 


Hie  development  of  the  explosion  can  be  viewed  in  terms  of  four  stages. 
These  are  (l)  the  spherical  free-field  expansion  prior  to  reflection  at  the 
outer  boundary;  (2)  multiple  distinct  reflections  off  both  boundaries,  during 
which  the  transition  to  one-dimensional  flow  takes  place;  (3)  approximately 
one  dimensional  expansion  after  reflecting  gas-dynamic  discontinuities  are  no 
longer  discernible;  (U)  a  gradual  transition  to  three-dimensional  flow  as  the 
finite  extent  of  the  array  comimnicates  itself  to  the  interior.  Stage  (l)  is 
not  contained  in  our  model;  since  the  weapon  mass  greatly  exceeds  the  mass  of 
air  contained  within  the  initial  explosion  radius,  the  nature  of  the  weapon 
and  its  constituents  ought  to  be  important  in  this  stage,  but  no  attempt  has 
been  made  to  improve  on  the  realism  of  the  1-kton  standard.  Likewise,  stage 
(U)  is  absent  from  the  model  (but  see  remarks  in  our  concluding  section 
below).  Instead,  we  attempted  to  model  stages  (2)  and  (3). 

The  calculation  was  performed  with  100  zones,  each  10.2  cm  in  width  in 
the  radial  direction,  and  100  zones  in  the  axial  direction,  90  of  which  were 
20. U  cm  in  length  and  the  final  10  of  which  geometrically  increased  by 
increments  of  11%  (Fig.  2a).  The  final  grid  size  was  10.2  m  by  22.1*  m.  An 
equivalent  18  Mton  surface  burst  was  achieved  by  inserting  values  of  density, 
energy,  and  momentum  scaled  from  the  1  kton  standard.  Since  values  of  the  1 
kton  standard  were  only  available  from  10  meters  (an  18  Mton  surface  burst 
scales  to  about  eight  meters  after  we  double  the  yield  to  account  for  surface 
reflection) ,  we  created  initial  profiles  by  increasing  the  density,  internal 
energy,  and  momentum  by  a  factor  of  1.1*.  As  a  result  the  initial  energy  was 
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3.1  x  1019  ergs  and  the  mass  was  8.2  x  106  g,  deposited  into  a  hemisphere  of 
10  m  radius.  Hie  effect  of  this  scaling  was  to  siirulate  the  initiation  of  an 
event  with  the  appropriate  yield  and  weapon  mass,  but  with  profiles  (as  a 
function  of  distance  from  the  burst  site)  which  were  not  correct  in  detail. 
During  the  early  phases  of  the  problem  they  must  give  rise  to  discrepant 
behavior,  but  at  later  times  this  presumably  becomes  insignificant. 

The  numerical  simulation  was  performed  using  the  NRL  FAST2D  code  which 
utilizes  the  latest  version  of  Flux  Corrected  Transport  (FCT).  The  code  was 
essentially  the  same  configuration  used  previously  for  lOU-ft  and  50-ft  1-kton 
HOB  calculations.  A  reflection  condition  was  imposed  on  the  left  and  bottom 
boundaries,  as  before,  and  on  the  right-hand  boundary  as  well.  In  addition  to 
the  diagnostics  employed  previously,  we  introduced  several  devised 
specifically  for  this  problem. 

Figure  2b  shows  plots  of  the  initial  profiles  of  the  mass  density, 
velocity,  pressure  and  energy  density.  Note  that  all  of  these  variables 
achieve  their  maxima  at  the  leading  shock.  To  convert  these  profiles  into 
their  counterparts  for  the  18  Mton  case,  it  is  necessary  to  dilate  by  a  factor 
of  (900  ft/10  m)1/3  =  27. 1*.  Likewise,  the  time  scale  has  to  be  expanded  by 
the  same  factor. 

The  reflection  of  this  initial  blast  wave  from  the  outer  radius  of  the 
can  and  the  subsequent  evolution  of  the  system  are  depicted  in  Fig.  3.  This 
shows  pressure  contours  and  velocity  vectors  plotted  at  intervals  of  100 
cycles  (approximately  30-50  us).  Although  the  time  intervals  vary  because  of 
the  changing  time  step,  these  plots  are  in  a  sense  equally  separated  in  terms 
of  the  "total  change"  in  the  system,  because  when  flow  and  sound  speeds  are 


large  the  timestep  decreases,  and  vice  versa.  Note  how  rapidly  the  reflected 

shock  propagates  to  the  left;  almost  invisible  at  cycle  201,  it  is  nearly 

halfway  across  the  mesh  at  cycle  301.  It  travels  upward  almost  as  rapidly, 

reaching  an  altitude  of  10  m  by  cycle  601.  By  cycle  801  it  appears  to  have 

overtaken  the  primary  shock  wave,  which  is  considerably  distorted  after  cycle 

1501.  On  the  right  hand  boundary  the  reflection,  which  is  initially  regular, 

gradually  converts  to  Mach  reflection.  This  development  is  of  course 

analogous  to  the  transition  to  Mach  reflection  observed  in  the  HOB  case3.  It 

is  hard  to  see  exactly  when  transition  occurs,  but  a  Mach  stem  is  clearly 

visible  in  the  pressure  plots  at  cycle  (1*0)  and  thereafter.  If  we  examine 

Carpenter’s  results  regarding  the  boundary  between  regular  and  Mach  reflection 
in  real  air  at  Mach  numbers  M  >  10,  we  see  that  the  angle  at  which  transition 

occurs  depends  on  M  only  weakly  and  is  approximately  equal  to  1*5°.  The  Mach 
stem  ought  thus  to  have  first  formed  around  cycle  701  to  801,  but  the  region 


of  reflection  is  not  sufficiently  well  resolved  to  show  it. 

In  Fig.  1*  another  view  of  the  same  time  sequence  is  shown.  These  surface 
pressure  plots  are  particularly  well  suited  for  showing  shock  and  rarefaction 
waves  propagating  through  the  sytem.  Since  the  pressures  are  interpolated 
onto  a  regular  mesh  before  plotting,  peak  values  are  lowered  by  as  much  as  a 
factor  of  two.  This  reduces  contrast  somewhat  and  makes  some  features  appear 
to  move  up  and  down  erratically  in  time,  but  the  overall  morphology  is  very 
clear.  We  see,  for  example,  how  the  shock  reflected  radially  from  the  origin 
appears  to  run  out  of  gas  at  about  cycle  901,  while  a  train  of  waves 
reverberates  down  from  the  top.  One  of  these  finally  reaches  the  radially- 
propagated  ground  shock  at  about  cycle  1801  and  begins  to  push  it  along. 
Throughout  this  time  the  pressures  in  the  lower  right  corner  of  the  system 
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remain  low.  Meanwhile,  the  shock  front  rising  toward  the  top  of  the  system, 
which  is  madly  churning  around  at  cycle  1201,  becomes  smoother  as  shocks 
propagate  back  down  from  it. 

Figure  5  shows  the  pressure  recorded  at  a  series  of  sensors  located 
across  the  bottom  of  the  can,  ranging  from  the  burst  center  (a)  to  the  outer 
edge  (t).  Note  that,  although  the  initial  profiles  have  a  sharp  (one-zone) 
discontinuity  at  the  leading  edge,  the  pressure  peak  at  sensor  l  takes  about 
TO  us  to  build  up  to  the  maximum.  This  is  about  a  hundred  times  as  long  as  it 
takes  the  shock  to  propagate  across  one  zone.  It  is,  however,  not  out  of 
line  with  the  time  required  for  the  air  behind  the  shock  (velocity  1.5  x  106 
cm/s)  to  cross  a  zone,  about  7  us.  The  density  (and  other  fluid  variables) 
can  build  up  to  their  post-reflection  values  when  the  air  from  (y+1)/(y-1) 

*  10  cells  has  been  confessed  into  the  cell  closest  to  the  wall. 

Using  Carpenter's4  reflection  factors  for  normally  incident  shocks,  we 
calculate  that  the  free-field  pressure,  initially  3.5  kbar,  should  reflect  up 
to  about  U0  kbar,  which  is  to  be  compared  with  the  value  of  25  kbar  we 
actually  obtained.  If  we  look  at  sensors  closer  to  the  burst  site,  we  see 
that  the  pressure  peaks  get  continuously  sharper,  finally  turning  into  an 
almost  discontinuous  spike  at  the  origin,  as  the  imploding  shock  attempts  to 
develop  into  a  singularity. 

The  information  in  these  pressure  histories  is  capsulized  in  Fig.  6, 
which  shows  the  maximum  pressure  recorded  during  the  course  of  the  calculation 
as  a  function  of  sensor  location.  Note  that  the  peak  at  r  =  0  is  nearly  as 
high  as  that  at  the  periphery,  in  contrast  with  results  reported  by  Fyatt5. 
This  may  be  attributable  to  the  lower  resolution  he  used,  or  to  the  greater 
numerical  dissipation  in  the  HULL  code.  Only  one  peak  was  detected  at  each 
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sensor  location*  Our  calculation  was  run  for  less  than  1  ms,  or  (scaled  to  18 
Mton)  rather  less  than  the  30  ras  lyatt  required  to  see  a  second  peak  at  the 
periphery  when  the  re-reflected  shock  impinges  there;  the  two  calculations 
are  thus  not  inconsistent  in  this  regard. 

The  approach  to  one-dimensional  motion  may  be  inferred  from  the  contour 
plots  of  Fig.  4,  where  it  is  seen  that  the  leading  edge  of  the  upward- 
propagating  blast  wave  becomes  more  and  more  horizontal,  and  from  the  velocity 
vector  plots  of  Fig.  3,  which  show  that  the  flow  becomes  increasingly  vertical 
and  that  velocities  well  behind  the  upper  front  tend  to  die  away.  A  more 
quantitative  comparison  is  possible  if  we  compute  the  average  pressure  on  the 
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P  = 


2*  /  p(r,0)rdr 
0 _ 

2*  /R  rdr 

n 


(3) 


as  a  function  of  time.  The  result  is  shown  as  the  solid  line  in  Fig.  7.  For 
comparison,  formula  (l)  has  been  added  (broken  line),  using  yield  W  =  1.4  kton 
and  area  A  =»  100*  ra2,  with  y  =  1.2.  the  agreement  between  the  two  is 
striking,  particularly  since  we  expect  to  observe  it  only  asymptotically.  It 
is  clear  that  the  one-dimensional  approximation  becomes  valid  very  early, 
probably  because  the  reverberating  shocks  behind  the  leading  front  of  the 
blast  propagate  so  rapidly  in  the  hot  medium. 


3.  CONCLUSIONS 


We  have  successfully  modelled  the  "bomb-in-a  can"  using  FAST2D,  a  general- 
purpose  code  which  required  no  special  modification  for  this  problem. 

Although  we  only  ran  the  calculation  out  to  an  altitude  of  about  20  m  (twice 
the  system  radius),  the  sliding  grid  (continuous  adaptive  rezone)  capability 
of  FAST2D  allows  us  to  continue  the  run  as  long  as  desired  by  stretching  the 
mesh  in  the  vertical  direction. 

Our  results  show  that  the  asymptotic  one-dimensional  state  is  approached 
very  rapidly,  apparently  because  the  mltiply  reflected  shocks  propagate  ruch 
faster  than  the  original  blast  wave.  Two  or  three  reflections  across  the 
radius  of  the  system  effectively  equilibrate  the  flow  and  relax  it  to  a  state 
of  vertical  expansion.  The  simulation  reveals  a  complex  pattern  of 
reverberations,  with  reflection  occurring  off  all  the  boundaries  and  regions 
where  conditions  are  highly  nonuniforra. 

A  number  of  modifications  are  required  to  make  these  results  applicable  to 
an  actual  tactical  situation.  Perhaps  the  most  inportant  is  the  inclusion  of 
dust  scoured  up  from  the  ground  and  entrained  in  the  wind  fields  following  the 
blast.  A  realistic  description  would  require  not  only  that  the  dust  mass  load 
the  air,  but  that  air  and  dust  be  permitted  to  exchange  momentum  and  energy. 
The  bottom  boundary  conditions  should  also  be  changed  to  mode*  the  energy  that 
goes  into  cratering  and  scouring. 

If  the  calculation  is  to  be  continued  to  late  times,  atmospheric 
stratification  should  be  included.  The  effects  of  water  vapor  and  turbulence 
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should  also  be  modelled.  Also,  venting  of  the  explosion  in  the  horizontal 
direction  (in  the  case  of  an  RV  exploding  near  the  edge  of  a  finite-eized 
array  of  shelters,  or  for  an  explosion  in  the  interior  of  the  array  which  has 
begun  to  feel  the  pressure  in  the  outer  explosions  drop  as  they  vent  outward) 
can  be  modelled  by  having  the  radius  R  increase  as  a  function  of  height.  Of 
course,  airblast  situations  in  which  the  explosion  originates  above  the  ground 
can  also  be  investigated  (in  which  case  Mach  stems  could  appear  on  the  ground 
as  well  as  the  sides  of  the  system). 

Although  some  of  these  extensions  require  nontrivial  code  modifications, 
there  is  no  reason  in  principal  why  the  present  results  could  not  be  augmented 
and  refined  greatly  by  additional  calculations. 
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Fig.  1  —  Schematic  of  unbounded  hexagonal  array  of  shelters  with  1800  ft 
i  separations  on  a  perfectly  reflecting  plane.  Shown  is  the  cross  section  of  the 

hexagonal  parallelpiped  formed  by  the  planes  bisecting  the  lines  connecting 
a  given  shelter  with  its  six  nearest  neighbors  (reflection  planes),  along  with 
the  coaxial  cylinder  having  the  same  cross  sectional  area. 
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Fig.  3a  —  Pressure  contours  and  velocity  vectors  calculated  using  FAST2D,  shown  at  intervals  of  100  cycles 
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Fig.  3f  —  Pressure  contours  and  velocity  vectors  calculated  using  FAST2D,  shown  at  intervals  of  100  cycles 
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Pressure  contours  and  velocity  vectors  calculated  using  FAST2D,  shown  at  intervals  of  100  cycles 
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Fig.  4a  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4b  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4c  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 


Fig.  4d  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar . 
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Fig.  4f  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 


Fig.  4g  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4h  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4k  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection 
The  pressure  scale  here  is  22  kbar. 
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Fig.  41  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4n  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4p  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4q  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection 
The  pressure  scale  here  is  22  kbar. 


Fig.  4r  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4s  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4t  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  4u  —  Pressure  contours  of  Fig.  3,  shown  in  orthographic  projection. 
The  pressure  scale  here  is  22  kbar. 
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Fig.  5a  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  0.0m 


Fig.  5b  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  0.035m 
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Fig.  5c  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  0.075m 
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Fig.  5d  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  1 .25m 
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Fig.  5e  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  2.0m 
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Fig.  5f  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  3.0m 
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Fig.  5g  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  4.5m 
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Fig.  5h  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  6.05m 
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Fig.  5k  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  9.5m 


Fig.  51  —  Pressure  histories  as  measured  by  sensors  located  at  radii  of  10.0m 
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Fig.  7  —  Average  pressure  on  floor  of  mesh  (solid  curve)  and  pressure  at  origin 
calculated  from  Sedov  one-dimensional  self-similar  solution  with  y  =  1.2 
(broken  curve)  as  functions  of  time 
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